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In this chapter we review the characterization of entanglement in Gaus- 
sian states of continuous variable systems. For two-mode Gaussian states, 
we discuss how their bipartite entanglement can be accurately quan- 
tified in terms of the global and local amounts of mixedness, and ef- 
ficiently estimated by direct measurements of the associated purities. 
For multimode Gaussian states endowed with local symmetry with re- 
spect to a given bipartition, we show how the multimode block entan- 
glement can be completely and reversibly localized onto a single pair 
of modes by local, unitary operations. We then analyze the distribution 
of entanglement among multiple parties in multimode Gaussian states. 
We introduce the continuous-variable tangle to quantify entanglement 
sharing in Gaussian states and we prove that it satisfies the Coffman- 
Kundu-Wootters monogamy inequality. Nevertheless, we show that pure, 
symmetric three-mode Gaussian states, at variance with their discrete- 
variable counterparts, allow a promiscuous sharing of quantum corre- 
lations, exhibiting both maximum tripartite residual entanglement and 
maximum couplewise entanglement between any pair of modes. Finally, 
we investigate the connection between multipartite entanglement and the 
optimal fidelity in a continuous- variable quantum teleportation network. 
We show how the fidelity can be maximized in terms of the best prepara- 
tion of the shared entangled resources and, viceversa, that this optimal 
fidelity provides a clearcut operational interpretation of several measures 
of bipartite and multipartite entanglement, including the entanglement 
of formation, the localizable entanglement, and the continuous-variable 
tangle. 
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1. Introduction 

One of the main challenges in fundamental quantum theory as well as 
in quantum information and computation sciences lies in the characteriza- 
tion and quantification of bipartite entanglement for mixed states, and in 
the definition and interpretation of multipartite entanglement both for pure 
states and in the presence of mixedness. While important insights have been 
gained on these issues in the context of qubit systems, a less satisfactory 
understanding has been achieved until recent times on higher-dimensional 
systems, as the structure of entangled states in Hilbert spaces of high di- 
mensionality exhibits a formidable degree of complexity. However, and quite 
remarkably, in infinite-dimensional Hilbert spaces of continuous-variable 
systems, ongoing and coordinated efforts by different research groups have 
led to important progresses in the understanding of the entanglement prop- 
erties of a restricted class of states, the so-called Gaussian states. These 
states, besides being of great importance both from a fundamental point 
of view and in practical applications, share peculiar features that make 
their structural properties amenable to accurate and detailed theoretical 
analysis. It is the aim of this chapter to review some of the most recent 
results on the characterization and quantification of bipartite and multi- 
partite entanglement in Gaussian states of continuous variable systems, 
their relationships with standard measures of purity and mixedness, and 
their operational interpretations in practical applications such as quantum 
communication, information transfer, and quantum teleportation. 

2. Gaussian States of Continuous Variable Systems 

We consider a continuous variable (CV) system consisting of N canon- 
ical bosonic modes, associated to an infinite-dimensional Hilbert space Ti. 
and described by the vector X — {xi,pi, . . . ,xn,Pn} of the field quadra- 
ture ("position" and "momentum") operators. The quadrature phase op- 
erators are connected to the annihilation a.i and creation aj operators of 
each mode, by the relations Xi = {cii + aj) and pi = (a^ — aj)/i. The 
canonical commutation relations for the X^'s can be expressed in ma- 
trix form: = 2ir2y, with the symplectic form Q = (Sf^iUJ and 
uj = 6ij-i - 6ij+i, i,j = 1, 2. 

Quantum states of paramount importance in CV systems are the so- 
called Gaussian states, i.e. states with Gaussian characteristic functions 
and quasi-probabihty distribution^ The interest in this special class of 
states (important examples include vacua, coherent, squeezed, thermal, and 
squeezed-thermal states of the electromagnetic field) stems from the feasi- 
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bility to produce and control them with hncar optical elements, and from 
the increasing number of efficient proposals and successful experimental im- 
plementations of CV quantum information and communication processes 
involving multimode Gaussian states (see Ref. |2l for recent reviews) . By 
definition, a Gaussian state is completely characterized by first and sec- 
ond moments of the canonical operators. When addressing physical prop- 
erties invariant under local unitary transformations, such as mixedness and 
entanglement, one can neglect first moments and completely characterize 
Gaussian states by the 2N x 2N real covariance matrix (CM) a, whose en- 
tries are — l/2{{Xi, Xj}) — {Xi){Xj). Throughout this chapter, cr will 
be used indifferently to indicate the CM of a Gaussian state or the state 
itself. A real, symmetric matrix cr must fulfill the Robertson-Schrodinger 
uncertainty relatioJ^l 



to be a bona fide CM of a physical state. Symplectic operations (i.e. be- 
longing to the group Sp(^2N,R) = {S e SL{2N,R) : S'^QS = Q}) acting 
by congruence on CMs in phase space, amount to unitary operations on 
density matrices in Hilbert space. In phase space, any A^-mode Gaussian 
state can be transformed by symplectic operations in its Williamson di- 
agonal forrri^J u, such that cr — S^uS, with v> — diag {^i, vi, . . . vn, J^n}- 
The set E = {vi} of the positive-defined eigenvalues of \iQ(T\ constitutes 
the symplectic spectrum of cr and its elements, the so-called symplectic 
eigenvalues, must fulfill the conditions i^i > I, following from Eq. ([T]) and 
ensuring positivity of the density matrix associated to cr. We remark that 
the full saturation of the uncertainty principle can only be achieved by 
pure iV-mode Gaussian states, for which I'i ~ 1 Vz = 1,...,N. Instead, 
those mixed states such that Vi<k = 1 and Viyk > 1, with I < k < N, 
partially saturate the uncertainty principle, with partial saturation becom- 
ing weaker with decreasing k. The symplectic eigenvalues Vi are determined 
by N symplectic invariants associated to the characteristic polynomial of 
the matrix |ir2cr|. Global invariants include the determinant Det cr ~ Yli 
and the quantity A = uf, which is the sum of the determinants of all 
the 2x2 submatrices of cr related to each modJ^. 

The degree of information about the preparation of a quantum state g 
can be characterized by its purity p = Tr , ranging from (completely 
mixed states) to 1 (pure states). For a Gaussian state with CM a one haa^ 




As for the entanglement, we recall that positivity of the CM's partial 
transpose (PPT}^ is a necessary and sufficient condition of separability 



(T + ifi > 0, 



(1) 




fi = l/V Det cr . 



(2) 
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for (M + A^)-mode bisymmetric Gaussian states (see Sec. 4) with respect 
to the M\N bipartition of the modes-, as weh as for (M + A^)-niode Gaus- 
sian states with fully degenerate symplectic spectruni^. In the special, but 
important case M = 1, PPT is a nece ssary and sufficient condition for sep- 
arability of all Gaussian stateJ^^^^. For a general Gaussian state of any 
Af |7V bipartition, the PPT criterion is replaced by another necessary and 
sufficient condition stating that a CM cr corresponds to a separable state 
if and only if there exists a pair of CMs cta and ctb, relative to the sub- 
systems A and B respectively, such that the following inequality holdJUl 
cr > (J A (B (tb- This criterion is not very useful in practice. Alternatively, 
one can introduce an operational criterion based on a nonlinear map , that 
is independent of (and strictly stronger than) the PPT conditiorM^. 

In phase space, partial transposition amounts to a mirror reflection of 
one quadrature in the reduced CM of one of the parties. If {i>i} is the 
symplectic spectrum of the partially transposed CM it, then a (l-l-A^)-mode 
(or bisymmetric (Af -I- A^)-mode) Gaussian state with CM cr is separable 
if and only if i>i > IV i. A proper measure of CV entanglement is the 
logarithmic negativit'^^ £V = log||^|li, where || • ||i denotes the trace 
norm, which constitutes an upper bound to the distillable entanglement of 
the state g. It can be computed in terms of the symplectic spectrum i>i of 



E^f = max j 0, -'V log i>i \ . (3) 
Ej^ quantifies the extent to which the PPT condition i>i > 1 is violated. 



3. Two Mode Gaussian States: Entanglement and 
Mixedness 

Two-mode Gaussian states represent the prototypical quantum states 
of CV systems, and constitute an ideal test-ground for the theoretical and 
experimental investigation of CV entanglementfl^l. Their CM can be written 
is the following block form 



a 7 



(4) 



where the three 2x2 matrices a, /3, 7 are, respectively, the CMs of the two 
reduced modes and the correlation matrix between them. It is well knowrP^ 
that for any two-mode CM cr there exists a local symplectic operation 
Si = Si® 5*2 which takes cr to its standard form (Jsf , characterized by 

a. = diag{a, a}, (3 = diag{6, 6}, 7 = diag{c+, c_} . (5) 
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States whose standard form fulfills a = b are said to be symmetric. Any 
pure state is symmetric and fulfills c+ — —c- = y/a^ — 1. The uncertainty 
principle Ineq. ([T]) can be recast as a constraint on the 5'p(4.R) invariants 
Deter and A(cr) = Deta + Det^ + 2Det7, yielding A(cr) < 1 + Deter. The 
standard form covariances a, b, c+, and c_ can be determined in terms of 
the two local symplectic invariants 

Ail = (Det a)-^/^ = 1/a , = (Det (3)-^^^ 1/6 , (6) 

which are the marginal purities of the reduced single-mode states, and of 
the two global symplectic invariants 

fi = (Det cr)-i/2 ^ _ _ c^)]-i/2 ^ A^a^+b'' + 2c+c_ , (7) 

where /i is the global purity of the state. Eqs. ([BllZl) can be inverted to 
provide the following physical parametrization of two-mode states in terms 
of the four independent parameters ^i, /i2, and 



Ml M2 4 



with Ezp = Y [A — (/^i =F /X2)^/(//^/i2)] — 4//i^. The uncertainty principle 
and the existence of the radicals appearing in Eq. ([8]) impose the following 
constraints on the four invariants in order to describe a physical state 

fJ-ifJ-2 < fJ- < —. r , (9) 

+ Imi - M2I 

- + ^^^i^<A < 1 + 1. (10) 

The physical meaning of these constraints, and the role of the extremal 
states {i.e. states whose invariants saturate the upper or lower bounds of 
Eqs. ([9l fT0)) ) in relation to the entanglement, will be investigated soon. 

In terms of symplectic invariants, partial transposition corresponds to 
flipping the sign of Det 7, so that A turns into A = A — 4 Det 7 = — A + 
2/^1 + 2//i2- The symplectic eigenvalues of the CM cr and of its partial 
transpose cr are promptly determined in terms of symplectic invariants 



21^1 = At \/A2 - 4/a*2 , 2i>^ = A T V - 4/^2 , (11) 

where in our naming convention < in general, and similarly for the 
i'lp. The PPT criterion yields a state cr separable if and only if > 1. 
Since Z/+ > 1 for all two-mode Gaussian states, the quantity i)^ also 
completely quantifies the entanglement, in fact the logarithmic negativ- 
ity Eq. Q is a monotonically decreasing and convex function of i)^, 
— max{0, — log I'-}. In the special instance of symmetric Gaussian 
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states, the entanglement of formation^ Ep is also computable^ but, be- 
ing again a decreasing function of v- , it provides the same characterization 
of entanglement and is thus fully equivalent to in this subcase. 

A first natural question that arises is whether there can exist two-mode 
Gaussian states of finite maximal entanglement at a given amount of mixed- 
ness of the global state. These states would be the analog of the maximally 
entangled mixed states (MEMS) that are known to exist for two-qubit 
systemJiSl, Unfortunately, it is easy to show that a similar question in 
the CV scenario is meaningless. Indeed, for any fixed, finite global purity 
H there exist infinitely many Gaussian states which are infinitely entan- 
gled. However, we can ask whether there exist maximally entangled states 
at fixed global and local purities. While this question does not yet have a 
satisfactory answer for two-qubit systems, in the CV scenario it turns out 
to be quite interesting and nontrivial. In this respect, a crucial observa- 
tion is that, at fixed fj,, /zi and ^2, the lowest symplectic eigenvalue of 
the partially transposed CM is a monotonically increasing function of the 
global invariant A. Due to the existence of exact a priori lower and upper 
bounds on A at fixed purities (see Ineq. [T0|) , this entails the existence of 
both maximally and minimally entangled Gaussian states. These classes of 
extremal states have been introduced in Ref. 1191 ^.nd completely character- 
ized (providing also schemes for their experimental production) in Ref. 1151 
where the relationship between entanglement and information has been ex- 
tended considering generalized entropic measures to quantify the degrees 
of mixedness. In particular, there exist maximally and minimally entangled 
states also at fixed global and local generalized Tsallis p-entropies"'^"^. In 
this short review chapter, we will discuss only the case in which the puri- 
ties (or, equivalently, the linear entropies) are used to measure the degree 
of mixedness of a quantum state. In this instance, the Gaussian maximally 
entangled mixed states (GMEMS) are two-mode squeezed thermal states, 
characterized by a fully degenerate symplectic spectrum; on the other hand, 
the Gaussian least entangled mixed states (GLEMS) are states of partial 
minimum uncertainty {i.e. with the lowest symplectic eigenvalue of their 
CM being equal to 1). Studying the separability of the extremal states (via 
the PPT criterion) , it is possible to classify the entanglement properties of 
all two-mode Gaussian states in the manifold spanned by the purities: 



MiA*2 < A* < 




^1+^2—^1^2 ' 



separable; 
coexistence; 

entangled. 



(12) 
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Fig. 1. Classification of the entanglement for two-mode Gaussian states in the space 
of marginal purities fii^2 and normalized global purity fj,/fiifi2. All physical states lie 
between the horizontal plane of product states /i = fii/i^, and the upper limiting surface 
representing GMEMMS. Separable states (dark grey area) and entangled states are well 
distinguished except for a narrow coexistence region (depicted in black) . In the entangled 
region the average logarithmic negativity (see text) grows from white to medium grey. 
The expressions of the boundaries between all these regions are collected in Eq. II12I I. 

In particular, apart from a narrow "coexistence" region where both 
separable and entangled Gaussian states can be found, the separability of 
two-mode states at given values of the purities is completely character- 
ized. For purities that saturate the upper bound in Ineq. ([9]), GMEMS 
and GLEMS coincide and we have a unique class of states whose entan- 
glement depends only on the marginal purities /ii,2- They are Gaussian 
maximally entangled states for fixed marginals (GMEMMS). The maximal 
entanglement of a Gaussian state decreases rapidly with increasing differ- 
ence of marginal purities, in analogy with finite-dimensional systemi2QI. 
For symmetric states (/ii = fi2) the upper bound of Ineq. ([9|) reduces to 
the trivial bound fj, < 1 and GMEMMS reduce to pure two-mode states. 
Knowledge of the global and marginal purities thus accurately characterizes 
the entanglement of two-mode Gaussian states, providing strong sufficient 
conditions and exact, analytical lower and upper bounds. As we will now 
show, marginal and global purities allow as well an accurate quantification 
of the entanglement. Outside the region of separability, GMEMS attain 
maximum logarithmic negativity -E^/max while, in the region of nonvan- 
ishing entanglement (see Eq. (I12p ). GLEMS acquire minimum logarithmic 
negativity i?A^min- Knowledge of the global purity, of the two local purities, 
and of the global invariant A [i.e., knowledge of the full covariance matrix) 
would allow for an exact quantification of the entanglement. However, we 
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Fig. 2. Maximal and minimal log- 
arithmic negativities as functions 
of the global and marginal puri- 
ties of symmetric two-mode Gaus- 
sian states. The darker (lighter) sur- 
face represents GMEMS (CLEMS). 
In this space, a generic two— mode 
mixed symmetric state is repre- 
sented by a dot lying inside the nar- 
row gap between the two extremal 
surfaces. 



will now show that an estimate based only on the knowledge of the exper- 
imentally measurable global and marginal purities turns out to be quite 
accurate. We can in fact quantify the entanglement of Gaussian states with 
given global and marginal purities by the average logarithmic negativity 
Ej^ = {Ej^/niax + -£'A/'min)/2 Wc Can then also define the relative error SE_\f 

on E^r as 5i?AA(Ml,2, = (-EA^max - i^A/-min)/(i^A/-max + E^f min) ■ It iS CaSy 

to see that this error decreases exponentially both with increasing global 
purity and decreasing marginal purities, i.e. with increasing entanglement, 
falling for instance below 5% for symmetric states {fii = /X2 = fii) and 
II > fj,i. The reliable quantification of quantum correlations in genuinely 
entangled two-mode Gaussian states is thus always assured by the experi- 
mental determination of the purities, except at most for a small set of states 
with very weak entanglement (states with Ej^ < 1). Moreover, the accu- 
racy is even greater in the general non-symmetric case /ii ^ fi2, because the 
maximal achievable entanglement decreases in such an instance. In Fig. [2l 
the surfaces of extremal logarithmic negativities are plotted versus /i^ and 
/i for symmetric states. In the case /i = 1 the upper and lower bounds co- 
incide, since for pure states the entanglement is completely quantified by 
the marginal purity. For mixed states this is not the case, but, as the plot 
shows, knowledge of the global and marginal purities strictly bounds the en- 
tanglement both from above and from below. This analysis shows that the 
average logarithmic negativity E^ is a reliable estimate of the logarithmic 
negativity E^, improving as the entanglement increases. We remark that 
the purities may be directly measured experimentally, without the need for 
a full tomographic reco nstruction of the whole CM, by exploiting quan t urn 
networks techniquei^ or single-photon detections without homodyninJ221 
Finally, it is worth remarking that most of the results presented here 
(including the sufficient conditions for entanglement based on knowledge of 
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the purities), being derived for CMs using the symplectic formahsm in phase 
space, retain their vaHdity for generic non Gaussian states of CV systems. 
For instance, any two-mode state with a CM equal to that of an entan- 

JOQI 

gled two-mode Gaussian state is entangled as well^. Our methods may 
thus serve to detect entanglement for a broader class of states in infinite- 
dimensional Hilbert spaces. The analysis briefly reviewed in this paragraph 
on the relationships between entanglement and mixedness, can be general- 
ized to multimode Gaussian states endowed with special symmetry under 
mode permutations, as we will show in the next section. 



4. Multimode Gaussian States: Unitarily Localizable 
Entanglement 

We will now consider Gaussian states of CV systems with an arbitrary 
number of modes, and briefly discuss the simplest instances in which the 
techniques introduced for two-mode Gaussian states can be generalized and 
turn out to be useful for the quantification and the scaling analysis of CV 
multimode entanglement. We introduce the notion of bisymmetric states, 
defined as those {M + A^)-mode Gaussian states, of a generic bipartition 
M\N, that are invariant under local mode permutations on the M-mode 
and A^-mode subsystems. The CM cr of a {M + A^)-mode bisymmetric 
Gaussian state results from a correlated combination of the fully symmetric 
blocks (Tq,m and ct^n : 

<r=(y ^ ) , (13) 

where cr„A/ (ct^n) describes a M-mode (A^-mode) reduced Gaussian state 
completely invariant under mode permutations, and F is a 2M x 2 A'' real 
matrix formed by identical 2x2 blocks 7. Clearly, F is responsible for 
the correlations existing between the M-mode and the A^-mode parties. 
The identity of the submatrices 7 is a consequence of the local invariance 
under mode exchange, internal to the A/-mode and A^-mode parties. A 
first observation is that the symplectic spectrum of the CM cr Eq. (fT3|) 
of a bisymmetric (M -I- A^)-mode Gaussian state includes two degenerate 
eigenvalues, with multiplicities M —1 and A'^ — 1. Such eigenvalues coincide, 
respectively, with the degenerate eigenvalue iy~ of the reduced CM ct^m 
and the degenerate eigenvalue of the reduced CM cr^jv, with the same 
respective multiplicities. Equipped with this result, one can provJSlthat a 
can be brought, by means of a local unitary operation, with respect to the 
M\N bipartition, to a tensor product of single-mode uncorrelated states 
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and of a two-mode Gaussian state with CM cr^''. Here we give an intuitive 
sketch of the proof (the detailed proof is given in Ref. W) . Let us focus on 
the iV-mode block ct^n. The matrices iflcrpN and iftcr possess a set of 
— 1 simultaneous eigenvectors, corresponding to the same (degenerate) 
eigenvalue. This fact suggests that the phase-space modes corresponding 
to such eigenvectors are the same for (t and for <t ^jv . Then, bringing by 
means of a local symplectic operation the CM ct^n in Williamson form, 
any {2N — 2) x {2N — 2) submatrix of a will be diagonalized because the 
normal modes are common to the global and local CMs. In other words, 
no correlations between the M-mode party with reduced CM ct^m and 
such modes will be left: all the correlations between the M-mode and A^- 
mode parties will be concentrated in the two conjugate quadratures of a 
single mode of the A^-mode block. Going through the same argument for 
the M-mode block with CM ct^a/ will prove the proposition and show that 
the whole entanglement between the two multimode blocks can always be 
concentrated in only two modes, one for each of the two multimode parties. 

While, as mentioned, the detailed proof of this result can be found in 
Ref. [8|(extending the findings obtained in Ref. |24]for the case M = 1), here 
we will focus on its relevant physical consequences, the main one being that 
the bipartite M x N entanglement of bisymmetric (M-l-iV)-mode Gaussian 
states is unitarily localizable, i.e., through local unitary operations, it can be 
fully concentrated on a single pair of modes, one belonging to party (block) 
M, the other belonging to party (block) TV. The notion of "unitarily local- 
izable entanglement" is different from that of "localizable entanglement" 
originally introduced by Verstraete, Popp, and Cirac for spin systemi^Sl. 
There, it was defined as the maximal entanglement conccntrable on two 
chosen spins through local measurements on all the other spins. Here, the 
local operations that concentrate all the multimode entanglement on two 
modes are unitary and involve the two chosen modes as well, as parts of 
the respective blocks. Furthermore, the unitarily localizable entanglement 
(when computable) is always stronger than the localizable entanglement. 
In fact, if we consider a generic bisymmetric multimode state of a M|A^ 
bipartition, with each of the two target modes owned respectively by one 
of the two parties (blocks), then the ensemble of optimal local measure- 
ments on the remaining ( "assisting" ) M + N ~ 2 modes belongs to the set 
of local operations and classical communication (LOCC) with respect to 
the considered bipartition. By definition the entanglement cannot increase 
under LOCC, which implies that the localized entanglement (a la Ver- 
straete, Popp, and Cirac) is always less or equal than the original M x N 
block entanglement. On the contrary, all of the same M x N original bi- 
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partite entanglement can be unitarily localized onto the two target modes, 
resulting in a reversible, of maximal efficiency, multimode/two-mode entan- 
glement switch. This fact can have a remarkable impact in the context of 
quantum repeaterJ^ for communications with continuous variables. The 
consequences of the unitary localizability are manifold. In particular, as al- 
ready previously mentioned, one can prove that the PPT (positivity of the 
partial transpose) criterion is a necessary and sufficient condition for the 
separability of {M + 7V)-mode bisymmetric Gaussian stateJ^. Therefore, 
the multimode block entanglement of bisymmetric (generally mixed) Gaus- 
sian states with CM cr, being equal to the bipartite entanglement of the 
equivalent two- mode localized state with CM cr'^'?, can be determined and 
quantified by the logarithmic negativity in the general instance and, for all 
multimode states whose two-mode equivalent Gaussian state is symmetric, 
the entanglement of formation between the Af-mode party and the iV-mode 
party can be computed exactly as well. 

For the sake of illustration, let us consider fully symmetric 27V-mode 
Gaussian states described by a 2N x 2N CM cr^2jv. These states are triv- 
ially bisymmetric under any bipartition of the modes, so that their block 
entanglement is always localizable by means of local symplectic operations. 
This class of states includes the pure, CV GHZ-type states (discussed in 
Refs.ElEll) 

that, in the limit of infinite squeezing, reduce to the simultane- 
ous eigenstates of the relative positions and the total momentum and co- 
incide with the proper Greenberger-Horne-Zeilingei^^ (GHZ) states of CV 
systemJ22l. The standard form CM (t^2n of this particular class of pure, 
symmetric multimode Gaussian states depends only on the local mixedness 
parameter b = which is the inverse of the purity of any single- mode 

reduced block, and it is proportional to the single-mode squeezing. Exploit- 
ing our previous analysis, we can compute the entanglement between a 
block of K modes and the remaining 2N — K modes for pure states (in this 
case the block entanglement is simply the Von Neumann entropy of any of 
the reduced blocks) and, remarkably, for mixed states as well. 

We can in fact consider a generic 27V-modc fully symmetric mixed state 
with CM cr^2N , obtained from a pure fully symmetric (2 -f Q)-niodc state 
by tracing out Q modes. For any Q and any dimension N of the block 
[K < N), and for any nonzero squeezing {i.e. for any b > 1) one has that 
the state exhibits genuine multipartite entanglement, as first remarked in 
Ref. 1271 for pure states: each K-mode party is entangled with the remaining 
{2N — K)-mode block. Furthermore, the genuine multipartite nature of the 
entanglement can be precisely quantified by observing that the logarithmic 
negativity between the K-mode and the remaining (2A^ — _ft')-mode block is 
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Fig. 3. Hierarchy of block en- 
tanglements of fully symmetric 
2Af-mode Gaussian states oi Kx 
{2N - K) bipartitions {N = 10) 
as a function of the single-mode 
mixedness b, for pure states 
(solid lines) and for mixed states 
obtained from (2A'^ + 4)-mode 
pure states by tracing out 4 
modes (dashed lines). 



an increasing function of the integer K < N , as shown in Fig. [3)The opti- 
mal sphtting of the modes, which yields the maximal, unitarily localizable 
entanglement, corresponds to K — N/2 if N is even, and K — {N — l)/2 
if N is odd. The multimode entanglement of mixed states remains finite 
also in the limit of infinite squeezing, while the multimode entanglement of 
pure states diverges with respect to any bipartition, as shown in Fig. [3l For 
a fixed amount of local mixedness, the scaling structure of the multimode 
entanglement with the number of modes can be analyzed as well, giving rise 
to an interesting r esultP. Let us consider, again for the sake of illustration, 
the class of fully symmetric 2A^-mode Gaussian states, but now at fixed 
single-mode purity. It is immediate to see that the entanglement between 
any two modes decreases with N, while the N\N entanglement increases 
(and diverges for pure states as ^ oo): the quantum correlations be- 
come distributed among all the modes. This is a clear signature of genuine 
multipartite entanglement and suggests a detailed analysis of its sharing 
properties, that will be discussed in the next section. The scaling struc- 
ture of multimode entanglement also elucidates the power of the unitary 
localizability as a strategy for entanglement purification, with its efficiency 
improving with increasing number of modes. Finally, let us remark that the 
local symplectic operations nee ded for the unitary localization can be im- 
plemented by only using passivJ^Sl ^jij active linear optical elements such 
as beam splitters, phase shifters and squeezers, and that the original mul- 
timode entanglement can be estimated by the knowledge of the gl obal a nd 
local purities of the equivalent, localized two-mode state (see Refs. ^2lIHl for 
a thorough discussion), along the lines presented in Section [3] above. 



5. Entanglement Sharing of Gaussian States 

Here we address the problem of entanglement sharing among multi- 
ple parties, investigating the structure of multipartite entang 

lementPED. 
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Our aim is to analyze the distribution of entanglement between different 
(partitions of) modes in CV systems. In Ref.i32iCoffman. Kundu and Woot- 
ters (CKW) proved for a three-qubit system ABC, and conjectured for A'' 
qubits (this conjecture has now been proven by Osborne and Verstraet j^^l) ^ 
that the entanglement between, say, qubit A and the remaining two-qubits 
partition (BC) is never smaller than the sum of the A|B and A|C bipartite 
entanglements in the reduced states. This statement quantifies the so-called 
monogamy of quantum entanglemenlf^S, in opposition to the classical cor- 
relations which can be freely shared. One would expect a similar inequality 
to hold for three-mode Gaussian states, namely 



where E is a proper measure of CV entanglement and the indices k} 
label the three modes. However, an immediate computation on symmet- 
ric states shows that Ineq. (ITil) can be violated for small values of the 
single-mode mixedness b using either the logarithmic negativity or the 
entanglement of formation Ep to quantify the bipartite entanglement. This 
is not a paradosfSH; rather, it implies that none of these two measures is 
the proper candidate for approaching the task of quantifying entanglement 
sharing in CV systems. This situation is reminiscent of the case of q ubit 
systems, for which the CKW inequality holds using the tangle r but 
fails if one chooses equivalent measures of bipartite entanglement such as 
the concurrence^ {i.e. the square root of the tangle) or the entanglement of 
formation itself. Related problems on inequivalent entanglement measures 
for the ordering of Gaussian states are discussed in Ref. [36l 

We then wish to define a new measure of CV entanglement able to cap- 
ture the entanglement distribution trade-off via the monogamy inequality 
(fT4| . A rigorous treatment of this problem is presented in Ref. [30l Here 
we briefly review the definition and main properties of the desired measure 
that quantifies entanglement sharing in CV systems. Because it can be re- 
garded as the continuous-variable analogue of the tangle, we will name it, 
in short, the contangle. 

For a pure state of a (1 A^)-mode CV system, we can formally 
define the contangle as 



Er{ip) is a proper measure of bipartite entanglement, being a convex, in- 
creasing function of the logarithmic negativity £V, which is equivalent 
to the entropy of entanglement in all pure states. For pure Gaussian 
states IV') with CM o-p, one has Er{crP) — log^(l/^i — \Jl/ ^\ — 1), where 



E^m _Ei\i -E'\^ >0, 



(14) 



i?,(V')=log2||g||i, g=\il:m. 



(15) 
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III = l/VDetcTi is the local purity of the reduced state of mode 1, de- 
scribed by a CM eri (considering 1 x bipartitions). Definition (|15p is 
extended to generic mixed states g of {N + l)-mode CV systems through 
the convex-roof formalism, namely: 

EriQ)= inf y^p^Erii^,), (16) 

where the infimum is taken over the decompositions of g in terms of pure 
states For infinite-dimensional Hilbert spaces the index i is contin- 

uous, the sum in Eq. (|16p is replaced by an integral, and the probabilities 
{pi} by a distribution n^ip). All multimode mixed Gaussian states cr ad- 
mit a decomposition in terms of an ensemble of pure Gaussian states. The 
infimum of the average contangle, taken over all pure Gaussian decomposi- 
tions only, defines the Gaussian contangle Gt, which is an upper bound to 
the true contangle i?r, and an entanglement monoton e unde r Gaussian lo- 
cal operations and classical communications The Gaussian 
contangle, similarly to the Gaussian entanglement of formation"^-', acquires 
the simple form Gr(cr) = infCTP<CT Er{cr^), where the infimum runs over all 
pure Gaussian states with CM cr^ < <j. 

Equip ped with these properties and definitions, one can prove sev- 
eral resultJ^. In particular, the general (multimode) monogamy inequality 
£;i,„|(zi...i™_ii™+i...iiv) _ Y^^^^ E'^^\" > is satisfied by aU pure three-mode 
and all pure symmetric A^-mode Gaussian states, using either Er or Gt 
to quantify bipartite entanglement, and by all the corresponding mixed 
states using Gr- Furthermore, there is numerical evidence supporting the 
conjecture that the general CKW inequality should hold for all nonsymmet- 
ric A^-mode Gaussian states as welljf The sharing constraint ^4]) leads to 
the definition of the residual contangle as a tripartite entanglement quanti- 
fier. For nonsymmctric three-mode Gaussian states the residual contangle 
is partition-dependent. In this respect, a proper quantification of tripartite 
entanglement is provided by the minimum residual contangle 



= min 



(17) 



where {i,j, k) denotes all the permutations of the three mode indexes. This 
definition ensures that Er''^'' is invariant under mode permutations and 



^The conjectured monogamy inequality for all (pure or mixed) Af-mode Gaussian states 
has been indeed proven by considering a slightly different version of the continuous- 
variable tangle, defined in terms of the (convex-roof extended) squared negativity instead 
of the squared logarithmic negativity [T. Hiroshima, G. Adesso and F. lUuminati, Phys. 
Rev. Lett. 98, 050503 (2007)]. 
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is thus a genuine three-way property of any three-mode Gaussian state. 
We can adopt an analogous definition for the minimum residual Gaussian 
contangle Gt"''*^- One finds that the latter is a proper measure of genuine 
tripartite CV entanglement, because it is an entanglement monotone under 
tripartite GLOCC for pure three-mode Gaussian stateJ^. 

Let us now analyze the sharing structure of multipartite CV entangle- 
ment, by taking the residual contangle as a measure of tripartite entan- 
glement. We pose the problem of identifying the three-mode analogues 
of the two inequivalent classes of fully inseparable three-qubit states. 



l^pw) = (l/\/3) [|001) + |010) + |100)]. These states are both pure and fully 
symmetric, but the GHZ state possesses maximal three-party tangle with 
no two-party quantum correlations, while the W state contains the maxi- 
mal two-party entanglement between any pair of qubits and its tripartite 
residual tangle is consequently zero. 

Surprisingly enough, in symmetric three-mode Gaussian states, if one 
aims at maximizing (at given single-mode squeezing b) either the two-mode 
contangle Er^ in any reduced state {i.e. aiming at the CV VF-like state), 
or the genuine tripartite contangle {i.e. aiming at the CV GHZ-like state), 
one finds the same, unique family of pure symmetric thr ee-m ode squeezed 
states. These states, previously named "GHZ-type" statei^S^ have been in- 
troduced for generic iV-mode CV systems in the previ ous S ection, where 
their multimode entanglement scaling has been studiecJSlHl. 

The peculiar 

nature of entanglement sharing in this class of states, now baptized CV 
GUZ/W states, is further confirmed by the following observation. If one 

i\(jk) 

requires maximization of the 1x2 bipartite contangle Et under the con- 
straint of separability of all two-mode reductions, one finds a class of sym- 
metric mixed states whose tripartite residual contangle is strictly smaller 
than the one of the GUZ/W states, at fixed local squeezin^SS Therefore, 
in symmetric three-mode Gaussian states, when there is no two-mode en- 
tanglement, the three-party one is not enhanced, but frustrated. 

These results, unveiling a major difference between discrete- variable and 
CV systems, establish the promiscuous structure of entanglement sharing 
in symmetric Gaussian states. Being associated with degrees of freedom 
with continuous spectra, states of CV systems need not saturate the CKW 
inequality to achieve maximum couplewise correlations. In fact, without vio- 
lating the monogamy inequality (jl4p . pure symmetric three-mode Gaussian 
states are maximally three-way entangled and, at the same time, maximally 
robust against the loss of one of the modes due, for instance, to decoher- 
ence, as demonstrated in full detail in Ref. [39, This fact may promote 
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these states, experimentally realizable with the current technologjSO!^ as 
candidates for reliable CV quantum communication. Exploiting a three- 
mode CV GKZ/W state as a quantum channel can ensure for instance 
a tripartite quantum information protocol like a teleportation network or 
quantum secret sharing; or a standard, highly entangled two-mode chan- 
nel, after a unitary (reversible) localization has been performed through a 
single beam splitter; or, as well, a two-party quantum protocol with better- 
than-classical efficiency, even if one of the modes is lost due to decoherence. 
We will next focus on a relevant applicative setting of CV multipartite en- 
tanglement, in which various of its properties discussed so far will come in 
a natural relation. 

6. Exploiting Multipartite Entanglement: Optimal Fidelity 
of Continuous Variable Teleportation 

In this section we analyze an interesting application of multipartite CV 
entanglement: a quantum teleportation-network protocol, involving N users 
who share a genuine A^-partite entangled Gaussian resource, completely 
symmetric under permutations of the modes. In th e standard multiuser 
protocol, proposed by Van Loock and BraunsteirJS]^ two parties are ran- 
domly chosen as sender (Alice) and receiver (Bob), but, in order to accom- 
plish teleportation of an unknown coherent state. Bob needs the results of 
N — 2 momentum detections performed by the other cooperating parties. 
A nonclassical teleportation fidelity (i.e. > J^'^'-) between any pair of 
parties is sufficient for the presence of genuine TV-partite entanglement in 
the shared resource, while in general the converse is false (see e.g. Fig. 1 
of Ref. |4T|) . The fidelity, which quantifies the success of a teleportation ex- 
periment, is defined a,s T = {tp'^^\g°'^*\tp™), where "in" and "out" denote 
the input and the output state. reaches unity only for a perfect state 
transfer, = while without entanglement in the resource, by 

purely classical communication, an average fidelity of J^ci = 1/2 is the best 
that can be achieved if the alphabet of input states includes all coherent 
states with even weighlP^l. This teleportation network has been recently 
demonstrated experimentalljW^ by exploiting three-mode squeezed Gaus- 
sian states'^'^, yielding a best fidelity oiT — 0.64±0.02, an index of genuine 
tripartite entanglement. Our aim is to determine the optimal multi-user 
teleportation fidelity, and to extract from it a quantitative information on 
the multipartite entanglement in the shared resources. By "optimal" here 
we mean maximization of the fidelity over all local single-mode unitary op- 
erations, at fixed amounts of noise and entanglement in the shared resource. 
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We consider realistically mixed iV-mode Gaussian resources, obtained by 
combining a mixed momentum-squeezed state (with squeezing parameter 
ri) and iV — 1 mixed position-squeezed states (with squeezing parameter 
r2 ^ ri and in principle a different noise factor) into an Af-sphtteiSH (a 
sequence of — 1 suitably tuned beam splitters). The resulting state is a 
completely symmetric mixed Gaussian state of a TV-mode CV system. For 
a given thermal noise in the individual modes (comprising the unavoidable 
experimental imperfections), all the states with equal average squeezing 
f = (ri + r2)/2 are equivalent up to local single-mode unitary operations 
and possess, by definition, the same amount of multipartite entanglement 
with respect to any partition. The teleportation efficiency, instead, depends 
separately on the different single-mode squeezings. We have then the free- 
dom of unbalancing the local squeezings ri and r2 without changing the 
total entanglement in the resource, in order to single out the optimal form 
of the resource state, which enables a teleportation network with maximal 
fidelity. This analysis is straightforward (see Ref.[44]for details), but it yields 
several surprising side results. In particular, one finds that the optimal form 
of the shared A^-mode symmetric Gaussian states, for iV > 2, is neither un- 
biased in the Xi and pi quadratures (like the states discussed in Ref. [45] for 
iV = 3), nor constructed by N equal squeezers (n = r2 = f). This latter 
case, which has been implemented experimentally^ for = 3, is clearly 
not optimal, yielding fidelities lower than 1/2 for iV > 30 and f falling in 
a certain interval^. According to the authors of Ref. SH the explanation 
of this paradoxical behavior should lie in the fact that their teleportation 
scheme might not be optimal. However, a closer analysis shows that the 
problem does not lie in the choice of the protocol, but rather in the choice 
of the resource states. If the shared iV-mode squeezed states are prepared, 
by local unitary operations, in the optimal form (described in detail in 
Ref. |44|) . the teleportation fidelity is guaranteed to be nonclassical 

(see Fig. H]) as soon as f > for any TV, in which case the considered class 
of pure states is genuinely multiparty entangled, as we have shown in the 
previous sections. In fact, one can sho^ that this nonclassical optimal 
fidelity is necessary and sufficient for the presence of multipartite entangle- 
ment in any multimode symmetric Gaussian state used as a shared resource 
for CV teleportation. These findings yield quite naturally a direct operative 
way to quantify multipartite entanglement in A'^-mode (mixed) symme tric 
Gaussian states, in terms of the so-called Entanglement of TeleportatioM^, 
defined as the normalized optimal fidelity 

Et = max {0, ~ T^) / {l - ] , (18) 
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Fig. 4. Plot of the optimal fi- 
delity for teleporting an arbi- 
trary coherent state from any 
sender to any receiver chosen 
from N (N = 2, . . . , 50) par- 
ties, exploiting A'^-party entan- 
gled, pure symmetric Gaussian 
states as resources. A nonclassi- 
cal fidelity J^'^'^ > 0.5 is always 
assured for any N, if the shared 
entangled resources are prepared 
in their optimal form. 



going from (separable states) to 1 (CV GHZ/VF state). Moreover, one 
finds that the optimal shared entanglement that allows for the maximal 
fidelity is exactly the CV counterpart of the localizable entanglement, orig- 
inally introduced for spin systems by Verstraete, Popp, and CiracP^. The 
CV localizable entanglement (not to be confused with the unitarily local- 
izable entanglement introduced in Section [?]) thus acquires a suggestive 
operational meaning in terms of tclcportation processes. In fact, the local- 
izable entanglement of formation (computed by finding the optimal set of 
local measurements — unitary transformations and nonunitary momentum 
detections — performed on the assisting modes to concentrate the high- 
est possible entanglement onto Alice and Bob pair of modes) is a mono- 
tonically increasing function of Et: E^"" = /[(I - Et) / {I + Et)], with 
f{x) = ii±f^ log _ (iz£)! log ii^. For TV = 2 (standard two-user 

teleportatiorP^ the state is already localized and Ep'^ = Ep, so that Et 
is equivalent to the entanglement of formation Ep of two-mode Gaussian 
states. Remarkably, for = 3, i.e. for three-mode pure Gaussian resource 
states, the residual contangle Er''^'' introduced in Section 5 (see Eq. p?)) ) 
turns out to be itself a monotonically increasing function of Et'. 



E^:>\k = log2 



2V2£'t - [Et 



{Et - 1)^Et{Et + ^) + 1 



l)^EfTT_l , 
2 ^ 



1 



Et{Et + '^) + 1 



(19) 



The quantity Et thus represents another equivalent quantification of gen- 
uine tripartite CV entanglement and provides the latter with an operational 
interpretation associated to the success of a three-party teleportation net- 
work. This suggests a possible experimental test of the promiscuous sharing 
of CV entanglement, consisting in the successful (with nonclassical optimal 
fidelity) implementation of both a three-user teleportation network exploit- 
ing pure symmetric Gaussian resources, and of two-user standard teleporta- 
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tion exploiting any reduced two-mode channel obtained discarding a mode 
from the original resource. 

Besides their theoretical aspects, the results reviewed in this section are 
of direct practical interest, as they answer the experimental need for the best 
preparation recipe of an entangled squeezed resource, in order to implement 
quantum teleportation and in general CV communication schemes with the 
highest possible efficiency. 

7. Conclusions and Outlook 

We have reviewed some recent results on the entanglement of Gaussian 
states of CV systems. For two-mode Gaussian states we have shown how 
bipartite entanglement can be qualified and quantified via the global and 
local degrees of purity. Suitable generalizations of the techniques introduced 
for two-mode Gaussian states allow to analyze various aspects of entangle- 
ment in multimode CV systems, and we have discussed recent findings on 
the scaling, localization, and sharing properties of multipartite entangle- 
ment in symmetric, bisymmctric, and generic multimode Gaussian states. 
Finally, we have shown that many of these properties acquire a clear and 
simple operational meaning in the context of CV quantum communication 
and teleportation networks. Generalizations and extensions of these results 
appear at hand, and we may expect further progress along these lines in 
the near future, both for Gaussian and non Gaussian states. 

A good portion of the material reported in this chapter originates from 
joint work with our friend and colleague Alessio Serafini, whom we warmly 
thank for the joy of collaborating together. It is as well a pleasure to ac- 
knowledge stimulating exchanges over the last two years with Nicolas Cerf, 
Ignacio Cirac, Jens Eisert, Jaromfr Fiurasek, Ole Kriiger, Gerd Leuchs, 
Klaus M0lmer, Tobias Osborne, Matteo Paris, Eugene Polzik, Gustavo 
Rigolin, Peter van Loock, Frank Verstraete, David Vitali, Reinhard Werner, 
Michael Wolf, and Bill Wootters. 
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